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ABSTRACT 



Aims. We attempt to compile a complete, mass-selected sample of galaxies with low specific star-formation rates, and compare their 
properties with theoretical model predictions. 

Methods. We use the f (24fim) / f (K) flux ratio and the SED fitting to the 0.35 - 8.0yum spectral distribution, to select quiescent 
galaxies from z - 0.4 to z — 4 in the GOODS-MUSIC sample. Our observational selection can be translated into thresholds in 
specific star-formation rate M/M,, which can be compared with theoretical predictions. 

Results. In the framework of the well-known global decline in quiescent galaxy fraction with redshift, we find that a non-negligible 
fraction ^ 15 - 20% of massive galaxies with low specific star-formation rate exists up to z — 4, including a tail of "red and dead" 
galaxies with M/M, < 

specific star-formation rates, but are unable to provide a global match to our data. 
Conclusions. 



10 "yr '. Theoretical models vary to a large extent in their predictions for the fraction of galaxies with low 
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1. Introduction 

Understanding the formation and evolution of early-type galax- 
ies is a major goal of present-day cosmology, as well as a funda- 
mental benchmark for "ab-initio" theoretical models of galaxy 
evolution. 

According to several independent lines of evidence, the pop- 
ulation of massive galaxies has undergone major evolution dur- 
ing the epoch corresponding to the redshift range 1 .5 - 3, where 
the gala xy stellar mass function evolve d significantly at high 
masses (Bert a et all 120071: iFontana et all 12006) (F06 in the fol- 
lowing), massive galax ies s ettled onto the Hubb l e sequ ence (e.g., 
lAbraham et al.l d2007l) and iFranceschini et alj d2006l)) and t he 
red sequence appears in high z clusters (e.g. JZirm et all 2008). 

The nature of the physical processes responsible for this 
rapid rise remains unclear. A large number of massive 
10 n M o ) actively star-forming galaxies is c learly in place at 
z - 2 dDaddi et al.ll2004tlPapovich et alj|2007l) . Within this pop- 
ulation, the m ore massive galaxie s tend to be the more actively 
star-forming (Dad di et al.ll200 7b). at variance with trends mea- 
sured in the local Universe. 

At the same time, galaxies with low levels of star-formation 
rates at z - 1.5-2 have be en detected by im aging surveys 
base d on color criteria ( e.g., iDaddi et aT] 120041) or SED fit- 
ting IGrazian et al.l l2007h. and by spectroscopic observations 
of red galaxy sam ples dCimatti et al1 l2004: Sara cco et al.ll2005t 
iKriek et aTl l2006f) . These results have motivated the inclusion 
of efficient methods for providing a rapid assembly of mas- 
sive galaxies at high z (such as starburst during interactions) as 
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well as quenching of the SFR, most notably via AGN feed back 
dMenci et al.ll2006l:lBower et alfeOOetlHopkins et al.ll2008l) . 

Unfortunately, a detailed validation of the prediction of the- 
oretical models has been hampered so far by the lack of a sta- 
tistically well defined sample of early-type galaxies at high red- 
shift, and by the difficulty in defining a common criterium to 
identify early-type galaxies in the data. It is difficult to iso- 
late passively evolving galaxies from the wider population of 
intrinsically red galaxies at high redshift, which include also 
a (probably larger) fraction of dust-enshrouded star-forming 
galaxies. The two classes are indeed indistinguishable when 
selected by means of a single color criterium, such as the 
"ERO " classification (R - K > 4) dDaddi et alj|2000l: iMcCarfhyl 
1 2004 or the "DRG" o ne {J - K > 2) dFranx et al.l 12003: 
van Dokk umet al.l 120031) . The corresponding SEDs, however, 
are not degenerate and can be distinguished with more com- 
plex criteria, even when spectroscopy is not feasible. Some of 
these criteria adopt more colors, such as t he / - J/ J - K method 
proposed by Pozzetti & Mannucci (2000) and spectroscopically 
val idated b y Cimatt iet alj d2003l) . or the BzK method proposed 
bv lDaddi etalj d2004l) . Other methods rely on the spectral fit- 
ting of the overall SED, either by making use of the resulting 
rest-fram e colors, as in the c ase of the U — V/V — I criteria pro- 
posed bv lWuvts et al.l d2007|). or directly using the output of the 
SEP fitting process dArnouts et al.ll2007l: ISalimbeni et al.ll2008l: 
IGrazian et al]|2007l) . 

We used the method of using SED fitting output in our analy- 
sis of the GOODS-S data set, to disentangle t he different contrib- 
utors to the mass density at high redshift (Gra zian et alj [2007) 
and to describe the evolution in the lumin osity function of red 
galaxies up to z — 3 (Salimbeni et al. 20Ql). 
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In this work, we use additional information available from 
observations in the mid-IR 24-fim band, to allow a more careful 
selection of high redshift passively evolving galaxies, with the 
specific aim of comparing their number density with theoretical 
expectations. 

We use a revise d version of our GOODS-MUSIC sample 
dGrazian et alj|2006h . a 15 band multicolor (U-to-24/im ) cat- 
alog extracted from an area of 143.2 arcminutes squared within 
the GOODS-South public survey. The key improvements com- 
pared to our previous work are a revised IRAC photometry 
and the addition of 24fim flux data for all object s in the cata- 
log. F ull details are given in a companion paper (Santi nTet al.l 
2009, S09 in the following). Another important difference is 
that we include all objects detected in the 4.5/im image, down to 
OT45 < 23.5. In the following, we adopt a mass-selected sample, 
obtained by applying a mass threshold at M, > 7 x 1O 1O M to 
our photometric catalog based on a combined selection K < 23.5 
or mus < 23.5. This photometric sample is complete at this mass 
limit to z — 4, also for dust-absorbed star-forming galaxies (see 
F06). 

Observed and rest-frame magnitudes are in the AB sys- 
tem, and we adopt the A-CDM concordance model (Hq = 
lOkmls/Mpc, Q. M = 0.3 and Q A = 0.7). 



2. Quiescent galaxies 

2.1. The role of mid-IR observations 

In principle, the observed mid-IR flux is a powerful tool for dis- 
tinguishing between the two classes of red galaxies at high red- 
shift. Dust-absorbed star-forming galaxies are expected to be 
bright in the mid-IR, where most of the UV-light absorbed by 
dust is re-emitted. In contrast, passively evolving galaxies are 
expected to be far dimmer, since ordinary stars have low emis- 
sion at IR wavelengths. To quantify this criterium, we use the 
ratio f(24fmt)/f(K) between the observed flux at 24/im and the 
Zf-band flux. We show the f(24/jm)/f(K) for all galaxies in our 
sample in Fig.Q] in four redshift bins. To show how this color can 
help us to differentiate between the galaxy types, we computed 
the same f(24fxm)/f(K) to z - 4 for a set of local templates in- 
cludin g early-type, spiral, and starburst galaxies ( Pollett a~et al.l 
|2007|) . All templates were fc-corrected at the different redshifts. 
As shown in Fig. \T\ the loci populated by the different galaxy 
types in the U — V versus f (24 jjm) / f(K) plane are separated 
well, allowing us to distinguish between actively star-forming 
galaxies and those with moderate-to-low star formation. 

In Fig.[T] we show the observed f (24 jjm) / f(K) for the com- 
plete M* > 7 x 10 10 M sample, both for individual galaxies 
as well as its general distribution. Upper limits were assumed 
to be equal to the l<x upper limits in the 24/vm photometry. It 
is clearly shown in Fig. Q]that the f(24fxm)/f(K) ratio of star- 
forming galaxies increases with redshift, and that a significant 
population of starburst exists at z > 1, with colors similar to 
those of local, rarer LIRGs. On the other hand, few bright galax- 
ies at low and intermediate redshifts, with low star-formation 
rates are detected at low levels in the 24/mi image. 

The most striking result is that f (24 jj.m)/f(K) distribution 
is clearly bimodal, its minimum corresponding to a gap be- 
tween the loci of active and passive galaxies. The observed 
minima in the distribution of the f(24fmi)/f(K) ratio occurs at 
f(24/jm)/f(K) = 1,2,5 and 6 in the four redshift bins adopted 
in Fig. [1] Two warnings are required about this bimodality. First, 
it is enhanced by the upper limits at faint 24fim levels, the true 
distribution probably being wider even at lower f(24fj.m)/ f(K). 
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Fig. 1. Ratio of the observed fluxes at 24/im to that in the K 
band as a function of the rest-frame U — V color, in four redshift 
bins, for the M, > 7 x 10 10 M o sample. Upper limits correspond 
to objects fainter than 20//Jy in the 24jum image. Galaxies with 
age /t > 6 (as discussed in Sect. 3) are shown in mag enta. Large 
blue c ircles are the "obscured AGN" candidates of iFiore et al.l 
(120081) . Closed areas represent the range o f f (24 ym) I f(K) ob - 
tained by redshifting a set of local samples (Poll etta et al.l l2007): 
from bottom to top, early-type galaxies (black line), spirals (blue 
line), starbursts (red line). 



More importantly, the bimodality disappears if a deeper cata- 
log is used (for instance, the z < 26 selected version of the 
GOODS-MUSIC catalog). The high mass-limit of our sample, 
indeed, preferentially selects extremely red galaxies, as shown 
by their significantly red U — V color, which are those for which 
the bimodality is clearly evident. This is a key feature of our ap- 
proach, since the f(24/jm)/f(K) ratio provides a natural way of 
distinguishing between active and quiescent galaxies within the 
population of high-redshift, red galaxies. 

2.2. The selection on rest-frame quantities 

The criterium that we adopted to identify quiescent galaxies can- 
not be directly applied to the output of current theoretical mod- 
els for galaxy formation and evolution, which typically do not 
predict the mid-IR flux. The comparison becomes more natural 
when we convert the data into the quantities immediately pro- 
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Fig. 2. Ratio between the observed fluxes at 24yum and in 
the K band as a function of the Specific Star-formation Rate 
(M/M*), for the sample with M* > 7 x 10 10 M G . Galaxies with 
M/M* < 1(T 12 yr~' have been arbitrarily set to M/M* - 1(T 12 
yr -1 . Upper limits refer to objects fainter than 20/iJy in the 
24/im im age. Large blue ci rcles are the "obscured AGN" can- 
didates of Fior e et all d2008l) . The solid vertical line corresponds 
to the inverse of the age of the Universe at the centre of the red- 
shift bin. The dashed vertical line shows the threshold on M/M* 
adopted to classify "red and dead" galaxies. The horizontal solid 
line refers to the minimum in the observed distribution of the 
f(2Ajj.m) I f(K) ratio in the same redshift range. 



vided by these models, i.e., stellar mass or star-formation rates. 
We use the stellar masses M, provided by the SED fitting, of ac- 
curacy described at length in F06, and briefly in Sect. 3. To de- 
rive star-form ation rates for galaxies detected at 24/itn, we apply 
the method of Pap ovich et al.l d2007l) in converting the 24/im + 
rest- frame UV luminosi ty into a total star-formation rate, using 
the lDale & Heloul d2002) models with a (lowering) correction at 
high mid-IR fluxes. For galaxies below 20/iJy at 24//m, we use 
the SFR derived from th e SED fitting. As we show in S09 (see 
also iDaddi et al.l 1200 7b). these two SFR estimators agree rela- 
tively well, especially at low SFR levels. 

In Fig. [2] we show the correlation between the 
f(24fim)/f(K) ratio and M/M* in our sample. For simplicity, 
we only plot the two redshift ranges that are more populated in 
our sample: the same trend holds at higher and lower redshifts. 
This correlation is largely expected given the relationship 
between M and the 24//m flux on the one hand, and between 
M* and the K magnitude on the other. The key point, however, 
is that the correlation is so tight that it allows us to translate 
the (observational) criterium based on the / '(24 fim)/f(K) ratio 
(horizontal line in Fig. |2]i into a (model-oriented) cut to the 




to- 12 10" 

M/M. (yr- 



Fig. 3. Probability distribution of the M/M* ratio in the sample 
of "red and dead" galaxies at 1.5 < z < 2.5, averaged over the 
whole sample. The thick curve shows the probability distribu- 
tion of the SED fitting to the 14 bands SED, from the U band 
to the 8//m band. The thin curve shows the probability distribu- 
tion removing models that over predict the SFR with respect to 
the upper limit provided by the constraint at 24yum. The dashed 
vertical line shows the limit M/M* < 10 _1I yr _1 that we adopt to 
define "red and dead" galaxies. 



estimated M/M* (vertical line in Fig. [2j. As shown in Fig. [2] the 
two samples in practice coincide. 

Based on this evidence, we use the specific star-formation 
rate M/M* to select quiescent galaxies, which allows direct com- 
parison of the data with theoretical models. Given that M/M* is 
dimensionally the inverse of a timescale, a natural threshold of 
M/M* is the inverse of the age of the Universe at the correspond- 
ing redshift (?E/(z)) _1 . We define galaxies with M/M* < (f(/(z)) _1 
as "quiescent" in the following. Such a name is motivated by the 
fact that - by definition - these galaxies have experienced a major 
episode of star-formation prior to the observations Q. 

3. "Red and dead" galaxies 

3.1. Basic definitions 

Galaxies defined as "quiescent" by the criterium presented in the 
previous Sect, may still have a measurable amount of ongoing 
SFR (e.g., M - 10M G /yr for a M* a 10 n M o galaxy). It is 
therefore interesting to estimate the fraction of galaxies with low 
or negligible levels of SFR. As we show in the following, the 
physical properties of these objects, that we label "red and dead", 
constitute a sterner test to theoretical models. 

To select galaxies at high redshift on the basis of low star- 
formation rates, we need to complement our f(24/jm) / f(K) ratio 
data with output from the fitting analysis of the optical-IR ob- 
served SED. The lack of detection in the 24/im image provides 
only an upper limit to the ongoing SFR (see for instance Fig. [2J, 
and its exact level can be estimated from the SED fitting only. 
We applied the SED fitting t echnique following th e recipe de- 
scribed in several papers (see iFontana et al.l (|2004|) and F06 for 
details). The U-to-8 um photometry wa s com pared with a grid 
of models from the Bruzual & Charlotl (l2003h (BC03) spectral 
synthesis code, characterized by exponentially declining star- 
formation histories of timescale r for a set of ages, metallicities, 
and dust extinctions. For comparison with our previous work, 
and most of the literature, we used the standard Salpeter IMF 



1 Indeed, if M. = {M) pas , x ty(z), where {M) pasl is the star-formation 
rate averaged over the whole age of the Universe at the corresponding 
z, the requirement M/M, < (%(z)) _1 implies M < {M) pasl . 
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and the BC03 models. We verified however that adoption of the 
most recent version of the code incor porating the treatment of 
the post-AGB stars dBruzual A||2007|) does not change our re- 
sults significantly (see Salimbeni et al. 2009 for a preliminary 
analysis). In particular, only 1 (out of 144) objects classified as 
"red and dead" would not be classified as such with the new ver- 
sion. 

To select "red and dead" galaxies, we use a threshold of 
M/M* < 10~ n y r 1 (vertical dashed line in Fig. |2j. To some ex- 
tent, this threshold is arbitrary, the important point being that 
we adopt the same cut in theoretical models, to ensure proper 
comparison. Ho wever, we note that the s ame threshold was also 
used previously (Brinchmann et al. 2004) to separate active and 
quiescent galaxies. We also note that it corresponds to a thresh- 
old age /t > 6 between the fitted age and the star-formation 
exponential timescale r. This follows from the fact that, assum- 
ing exponentially declining star-formation histories, the specific 
star-formation rate can be computed analytically to be M/M* = 
(r(e (fl s f/T) - I))" 1 , yielding M/M* lO -11 ^ -1 for age/r =* 6, 
for small values of r. This coincidence allows us to compare 
the pure SED fitting approach with the additional analysis pro- 
yided by the 24^m data. A threshold to age I r was applied before 
dArnouts et alj|2007t ISahmbeni et"aT1l2008l: iGrazian et ai]|2007l) 
to broadly identify passively evolving galaxies, although with 
a less conservative value of age/r > 4. In Fig. [T] we plot the 
location in the U — V vs. f(24fim)/f(K) plane of all galaxies 
with age/r > 6 (magenta points). The classification scheme 
based purely on the optical-near IR SED fitting agrees well 
with the /(24/ira) / f{K) criterium, providing an important con- 
sistency check. The main exception are a number of "Compton 
thick AGN " candi dates at 1.5 < z < 2.5, selected as described in 
iFiore et al] d2008f) . where the 24/jm emission was attributed to 
AGN activity. These objects are discussed further in Sect. 3.4. 

We select "red and dead" galaxies by requiring that M/M* < 
10 _11 yr _1 , not only that age/r > 6. This is because when es- 
timating SFR, we assign the SED-derived values only to ob- 
jects undetected in the mid-IR, i.e., those with 24/um flux below 
20/Jy. For this reason, the few objects fitted with age/r > 6 but 
with large f(24fj.m)/f(K) are not included in our sample of "red 
and dead" galaxies. We now discuss the accuracy of this method 
for the most interesting subsample, i.e., "red and dead" galaxies 
atz > 1.5. 

3.2. "Red and dead" galaxies at 1.5 < z < 2.5 

A sizeable fraction of the sample with M/M* < 10 _I1 yr _I 
have redshifts in the range 1.5 < z < 2.5, where we detect 
32 galaxies. The reliability of their M/M* measure can be as- 
sessed ijjshig_the_error analysis already widely adopted in similar 
cases dPapovichetal .1120011 F06). Briefly, the full synthetic li- 
brary used to determine the best-fit model spectrum is compared 
with the observed SED of each galaxy. For each spectral model 
(i.e., for each combination of the free parameters age, t, Z and 
E(B - V)) the probability P of the resulting x 2 is computed and 
retained, along with the associated M and M,. The ensemble of 
M/M* values corresponding to probabilities above a given mini- 
mum threshold defines the range of acceptable values for M/M*. 
Since most of our "red and dead" galaxies have only a photomet- 
ric redshift, our analysis for these objects has been completed by 
allowing the redshift parameter to be free. In Fig. [3] we show the 
resulting distribution of the probability P, averaged over the en- 
tire sample of "red and dead" galaxies at 1.5 < z < 2.5. The tail 
at large values of M/M* is due to 2 objects with a nearly equal 
probability of being fitted accurately by a dusty starburst model 



solution. The correspondingly large values of M that would be 
required would easily overpredict the observed flux (which is an 
upper limit) at 24/mi. Taking into account this additional con- 
straint, we display in Fig.[3]the distribution of P after removing 
the models that overpredict the SFR with respect to the value 
provided by the upper limit in the emission at 24/rni. It is clear 
from Fig. [3] that the SED fitting for these galaxies is well con- 
strained at levels of M/M* < 10 _11 yr _1 , since the probability of 
having M/M* has only a small tail above 10 n yr 1 . 

This result is unsurprising. Since the galaxies in our sam- 
ple are by selection massive (M« > 7 x 1O 1O M ), the thresh- 
old M/M* = 10~ n yr~' is equivalent to small but measur- 
able amounts of SFR, of approximately lM yr~'. These lev- 
els correspond to detectable fluxes in the deep B-z optical im- 
ages of GOODS (the magnitude computed from a galaxy with 
^ lM yr _1 , for small values of E(B - V), is in the range 26-28 
mags) and can then be measured by the full SED fitting. 

We can also check the small subset of our "red and dead" 
sample with spectroscopic observations in this redshift range. 
Using both the K20 and the public GOODS data set, we found 
6 objects that have a spectroscopic redshift. Four of these ob- 
jects are classified as early spectral type, with no detectable [Oil] 
emission line. Two other objects have both absorption features, 
typical of evolved systems, and a weak [Oil] line. A weak [Oil] 
line is not nece ssarily an indicator of ongoing star -formation 
dYan et al.l2006l) . Adopting the standard conversion of lKennicuttl 
d 1998b . we derive star formation rates of the order of 2-3 M Q yr~ 1 , 
consistent with the SED fitting estimates, confirming that these 
objects have M/M* < 10 _11 yr _1 . We conclude that our analy- 
sis is sufficiently robust at 1.5 < z < 2.5, and allows to build a 
mass-selected sample of "red and dead" galaxies. 

3.3. "Red and dead" galaxies atz> 2.5 

As we move to higher redshifts, 2.5 < z < 4, our analysis be- 
comes more uncertain. First, the constraints on the 24/im emis- 
sion are weaker or non-existent, since this band is shifted out- 
side the range dominated by dust emission. Objects also become 
even fainter and redder, and in general lack any spectroscopic 
confirmation. In this redshift range, we detect 12 "red and dead" 
candidates, 9 located at 2.5 < z p i, t < 3. We carefully exam- 
ined the data for each object to assess the detection reliability. 
First, we verified that the photometry was not contaminated by 
nearby companions, and that the overall SED was smooth and 
not biased obviously by photometric fakes. All galaxies selected 
in this redshift range are very red with z - B > 4, which is redder 
than typical EROs at lower redshifts. Their optical-IR SED is 
dominated by a break between the K and the IRAC bands, and 
a change of slope that is a signature of an evolved galaxy pop- 
ulation (star-forming, dusty objects exhibit a more featureless 
slope). The redshift probability distribution is broad, but typi- 
cally implies that z p i w t > 1.5-2 with a relatively large spread 
Az p hot - 0.5. 

Quantitatively, the error analysis described above indicated 
that 5 of the 9 objects at 2.5 < z p hot < 3 are constrained to have 
M/M* < 10" u y r ~\ even for such a broad redshift range. An 
example of these objects is shown in the upper panel of Fig. [4] 
The four other objects have a wider distribution of M/M (Fig. [4] 
middle panel), as well as redshift. 

In a similar way, the physical properties and the nature of the 
three objects detected at z > 3 are weakly constrained. They have 
a clear minimum in the^ 2 distribution around z p hot - 3.5-4, but 
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Fig. 4. Examples of "red and dead" galaxies at z > 2.5. For each 
object, from left to right: the observed flux in the GOODS band 
and the best-fit SED; the redshift probability function; the M/M* 
probability distribution. From top to bottom, the three galaxies 
shown represents three different categories: a robust candidate at 
2.5 < z < 3; a less secure candidate in the same redshift range; a 
typical example of the candidates at z > 3. 



have a tail to z — 2. The distribution of acceptable M/M* values 
also extends significantly beyond 10 _1I yr _1 . 

We conclude that at least 55% of our "red and dead" candi- 
dates at 2.5 < Zphoi < 3 are robust candidates, while the remain- 
ing 45% and the three candidates at z > 3 are more uncertain. 
It would be natural to convert this finding into an upper limit to 
the true number density of "red and dead" galaxies. However, 
we note that a comparable number of galaxies are inferred to 
have M/M* slightly above 10 _1I yr _1 and are not included in our 
sample, but have a range of acceptable model fits extending well 
below M/M* = 10 _1I yr _1 . We conclude that at z > 2.5 the es- 
timate of M/M* becomes unreliable, due mainly to the limita- 
tions in the depth of existing observations, and that the resulting 
statistical analysis of "red and dead" galaxies must be treated 
with caution. As we show, the statistical error (including both 
Poisson and cosmic-variance components) is so large for these 
small samples that a more rigorous treatment of these uncertain- 
ties is probably unnecessary. 

3.4. Hosts of Compton thick AGNs 

A potential source of uncertainties in many statistical analysis 
of high redshift galaxies is contamination by AGNs. As we de- 
scribe in S09, we removed from our catalog both spectroscopi- 
cally confirmed AGNs and X-ray sources with an optically dom- 
inant point-like source. We also identified galaxies that likely 
harbor Compton thick AGNs using the criterium defined by 
iFiore et a"D (120081) . which is based on the detection of a mid-IR 
excess in v ery red galaxies. T his criterium is similar to the one 
adopted bv lDaddi et al.l d2007al) but also includes objects that are 
much redder in the optical-IR range than the BzK-selected sam- 



ple of lDaddi etall (l2007al) . A full discussion of the SED of these 
objects will be presented elsewhere: for the moment, we note 
that 9 of these "Compton thick AGN" candidates at 1 .5 < z < 2.5 
are classified as "red and dead" galaxies, as shown in Fig. Q] 
Unfortunately, the nature of these objects remains elusive. Most 
of these objects are among the reddest in our sample, as indi- 
cated by their extreme rest-frame U — V. They are extremely 
faint or even undetected in the deep z-band ACS images, and 
are included in our sample only because of their detection in ei- 
ther K or 4.5^m images. As a result, the probability distributions 
of their photometric redshifts are often broad, although they are 
typically constrained to be at z p hot > 2. Only three (of nine) 
galaxies have the distribution of M/M* constrained to be below 
10" n y r_1 - In the following, we do not rely on this classifica- 
tion, and we therefore compute the fraction of "red and dead" 
galaxies at 1.5 < z < 2.5 in three different ways (labelled a), 
b), c) in the caption of Fig. |5), to explore all possible cases. 
First (option a)), we remove all "Compton thick AGN" candi- 
dates from our sample, irrespective of their SED classification. 
Alternatively (option b)), we assume that none of them is a "red 
and dead" galaxy, which provides a lower limit to their fraction. 
Finally (option c)), we include them in our statistical analysis, 
assigning their formal SED classification. As we show in the 
following Sect., the fraction computed with these different as- 
sumptions changes its value by ^20%, still providing interesting 
constraints on theoretical models. 

4. Results and discussion 

We finally compute the fraction of "quiescent" and "red and 
dead" galaxies in different redshift ranges, following the def- 
initions provided in the last two Sections. These fractions are 
shown in Fig|5]as a function of redshift. The different numbers 
of "red and dead" galaxies at 1.5 < z < 2.5 reflect the differ- 
ent accounting methods of hosts of Compton thick AGNs. Error 
bars were computed by summing (in quadrature) the Poisson 
and cosmic variance error. The latter was computed by measur- 
ing the relative variance with in 200 samples bootstrappe d from 
the Millennium Simulation dKitzbichler & White! |2007|) . using 
an area as large as GOODS-S and applying the same selection 
criteria. 

We recall that these fractions are computed with respect to 
the total number of galaxies with M* > 7 x 10 10 M o in our sam- 
ple. Since we are interested in the ratio between galaxy classes, 
the impact of the exact choice of this threshold is measurable but 
not dramatic. We verified that by either increasing or decreasing 
the mass limit cut by a factor of 2, the fraction of both "qui- 
escent" and "red and dead" galaxies in different redshift ranges 
changes by about 0.1, in the data (where applicable) as well as 
in the theoretical models. This factor of 2 variation may also be 
produced by adopting a different IMF, and is similar to the un- 
certainty in the stellar mass estimate. 

Our analysis confirms the cosmological decrease in the num- 
ber density of massive early-type galaxies at high redshifts: qui- 
escent galaxies dominate the population of massive galaxies at 
z « 1, and become progressively less common at higher z. 
However, we note that a significant fraction of galaxies with low 
levels of SFR is in place even at the highest redshifts sampled 
here (z - 3.5) with a fraction of about 10-15% at z > 2. 

By definition, "quiescent", massive galaxies assembled most 
of their stellar mass in previous epochs, implying that they expe- 
rienced an active starburst phase or important merging process at 
higher redshifts. The large fraction of "quiescent" galaxies that 
we observe at z - 2 implies that these processes must have been 
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Fig. 5. Lower Fraction of "red and dead" galaxies (defined by 
M/M, < lO^'yi" 1 ) as a function of redshift, in the M t > 
7 x 1O 1O M mass-selected sample. Points represent the observed 
values. Filled, open and starred points (slightly offset for clarity) 
refer to the three different strategies a,b,c discussed in the text to 
account for obscured AGNs. Error bars include Poisson and cos- 
mic variance errors. Lines refer to t he prediction s of t heoretical 
models, as describe d by the legend: Menci et alj (1200 6) (M06), 
Fontanot et al] d2007l) ( Km. iKitzbichler & White! d2007t) (K07). 
Nagamine et al. (2006), (N06); Upper Fraction of "quiescent" 
galaxies as a function of redshift, defined by M/M, < (?r/(z)) _1 , 
in the same mass-selected sample. Points and lines as in the 
lower panel. 



frequent at very high redshifts. The upper panel of Fig. [5] can 
also be interpreted in terms of "active" galaxies, i.e., the com- 
plementary galaxy fraction with M/M* > (£r/(z)) _1 . According 
to the duty-cycle argument presented in S09, these actively star 
forming galaxies experienced a major episode of star-formation, 
potentially building up a substantial fraction of their stellar mass 
in this episode. Their fraction increases with redshift, constitut- 
ing more than 50% of the massive galaxy population at z > 2. 

Finally, the sizeable number of "red and dead" galaxies is 
already in place at z > 1.5, implying that the star-formation 
episodes must be quenched either by efficient feedback mecha- 
nism and/or by the stochastic nature of the hierarchical merging 
process. 

It is interesting to determine whether theoretical models 
agree with these observational results. In Fig. [5] we plot the pre- 
dictions of several models, applying the same criteria based on 
the M/M, values. W e co nsider purely semi-ana lytical models 
dMenci et alj d2006h . M06. iFontanot et all d2007l) . F07), a semi- 
analytic al rendition of the Millenn ium N-body dark matter sim- 
ulation dKitzbichler & White! d2007b. K07) and purely hydrody- 
namical simulations (Nagami ne et alJ d2006). N06). The final are 
presented for three different timescales r of the star-formation 
rate (ranging from 2 x 10 7 yrs to 2 x 10 8 yrs), and represented 
with a shaded area. All these models agree in predicting a grad- 



ual decline with redshift in the fraction of galaxies with a low 
SFR. 

As far as the "quiescent" fraction is concerned, we note that 
most models agree quantitatively in their predictions at all red- 
shifts. There is a general tendency to slightly overpredict the 
fraction of "quiescent" galaxies, although the relatively large 
error bars of our sample prevent firm conclusions from being 
drawn . As we show in more detail in S09 (see also iDaddi et alj 
2007b), this is the result of an overall underestimate of the me- 
dian SSFR of massive galaxies, which increases the number of 
mildly star-forming galaxies that we detect within our selection 
criteria M/M* < (tuiz))" 1 - The F07 models provides a notice- 
able exception to this process. A main success of this model, 
indeed, is its capability of reproducing the S cuba counts and 
the high associated SFR (Fontan ot et al J 120071) : it is unsurpris- 
ing that it also predicts a large fraction of active galaxies, and 
hence a smaller fraction of "quiescent" ones. 

Conversely, the predicted fraction of "red and dead" galaxies 
varies significantly at all redshifts. This reveals that the predicted 
fraction of galaxies with very low levels of SFR is a particularly 
sensitive quantity, and provides a powerful way of highlighting 
the differences between the models. Some models (M06, F07) 
underpredict the fraction of "red and dead" galaxies at all red- 
shifts, and in particular predict virtually no object at z > 2, in 
contrast to what observed. The Millennium-based model agrees 
with the observed quantities, while the hydro model appears to 
overpredict them. 

It is beyond the scope of the present paper to discuss the ori- 
gins of these differences. They are difficult to ascertain, because 
of the complex interplay between all the physical processes in- 
volved in these models, the different physical process imple- 
mented - most notably those related to AGN feedback - and their 
different technical implementations. The failure of most models 
to reproduce simultaneously the fraction of "quiescent" and "red 
and dead" massive galaxies in the early Universe probably im- 
plies that the balance between the amount of cool gas and the 
star-formation efficiency on the one side, and the different feed- 
back mechanisms on the other, is still poorly understood. 
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